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ABSTRACT 


An  improved  computerized  technique  has  been  used  to 
analyze  nonstationary  cavitation  noise  from  an  oscil¬ 
lating  hydrofoil.  Combined  analog  and  digital  signal 
processing  procedures  were  developed  that  use  ensem¬ 
ble  averaging  techniques  to  quantify  the  average 
sound  pressure  variance  as  a  function  of  the  foil's 
angle  of  attack.  These  averages  provide  a  quanti¬ 
tative  basis  for  comparing  cavitation  noise  produced 
by  this  hydrofoil  under  a  finite  spectrum  of  specific 
test  conditions. 


INTRODUCTION 


Cavitation  is  a  destructive  phenomenon  associated  with  operating 
propellers  and  hydrofoils  when  the  local  pressure  on  their  lifting 
surfaces  is  less  than  the  vapor  pressure  of  the  water  —  the  water 
boils.  Cavitation  produces  noise,  metal  erosion,  loss  of  lift  or  thrust, 
and  vibration  in  surrounding  structures.  A  considerable  effort  is  underway 
in  the  Navy  to  eliminate  cavitation,  reduce  cavitation  to  acceptable 
levels,  or  to  counter  the  effects  of  cavitation.  For  example,  experiments 
are  performed  to  study  the  influence  of  the  static  and  time-varying 
pressure  fields  that  produce  cavitation.  Knowledge  obtained  from 
these  experiments  is  used  to  reduce  cavitation  through  enlightened 
design  of  propellers  and  hydrofoils,  and  to  predict  their  cavitation 
performance . 

Experiments  to  investigate  features  of  cavitation  noise  contribute 
useful  information  to  this  study.  One  such  experiment  with  an  oscillating 
hydrofoil  was  performed  in  the  David  W.  Taylor  Naval  Ship  Research 
and  Development  Center's  (DTNSRDC)  36  inch  Variable  Pressure  Water 
Tunnel*.^  During  this  experiment,  two  types  of  tests  were  conducted: 
static,  in  which  the  foil  was  fixed  at  a  given  angle  of  attack; 
and  dynamic,  in  which  the  foil's  angle  of  attack  was  sinusoidally 
varied  about  an  initial  offset.  Several  runs  of  each  type  of  test  were 


*A  complete  listing  of  references  is  given  on  page  27. 


made  at  different  water  speeds  and  tunnel  pressures.  Cavitation  noise 
data  were  recorded  and  photographs  of  the  cavitation  bubble  pattern 
were  taken . 


Cavitation  noise  was  statistically  quantified  for  each  of  three 
conditions: (1 )  the  sound  pressure  variance  at  the  instant  a  photograph 
was  taken,  (2)  the  average  variance  during  a  run,  and  (3)  the  average 
variance  as  a  function  of  the  hydrofoil's  angle  of  attack  during  dynamic 
runs.  Procedures  to  obtain  these  sound  pressure  measurements  for 
the  first  two  conditions  were  relatively  straightforward.  But  the 
third  was  not,  because  the  cavitation  noise  during  the  dynamic  runs 
was  modulated  by  the  foil's  motion  and  was  thus  statistically 
nonstationary. 

Typically  this  nonstationary  data  would  be  reduced  by  manually 
analyzing  oscillograph  records  of  the  modulated  noise.  This  experiment 
had  almost  100  dynamic  runs,  and  such  manual  analysis  would  have  consumed 
a  considerable  man-day  effort  and  period  of  time.  And  the  results 
from  manual  analysis  are  observer-dependent;  the  eye  and  mind  of  the 
individual  must  construct  the  image  of  the  "typical1'  noise  cycle. 

Accordingly,  special  analog  and  digital  signal  processing  procedures 
were  developed  to  automate  the  data  ^reduction  of  this  modulated  data 
and  to  obtain  observer-independent  quantitative  results.  Especially 
significant  was  the  ensemble-averaging  procedure  to  obtain  average 
sound  pressure  variance  as  a  function  of  the  foil's  angle  of  attack. 

This  report  includes  a  general  description  of  the  data  collected 
during  this  water  tunnel  experiment  as  well  as  a  detailed  description 
of  the  procedures  used  to  quantify  the  data.  The  first  section.  Test 
Data,  provides  an  account  of  how  the  experimental  data  were  acquired. 
Also,  in  this  section,  the  main  characteristics  of  the  hydrophone 
signal  (cavitation  noise  plus  contamination  noise)  are  related  to 
the  signal  processing  procedures  necessary  to  quantify  this  data. 

The  second  section.  Data  Reduction,  gives  a  detailed  description  of 
the  signal  processing  procedures  used  to  statistically  quantify  the 
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experimental  data*  First,  a  general  description  of  these  procedures 
is  presented;  then  the  procedures  are  detailed  in  two  parts:  Analog 
Signal  Processing  and  Digital  Signal  Processing.  The  last  section. 
Conclusions,  summarizes  the  quality  of  the  procedures  and  presents 
recommendations.  The  appendix  presents  a  sample  output  from  the  computer 
processing. 


TEST  DATA 


DATA  ACQUISITION 

The  experiment  took  place  in  DTNSRDC's  36  inch  Variable  Pressure  Water 
Tunnel.  A  hydrofoil  section  placed  in  the  flow  was  subjected  to  various 
water  speeds  and  tunnel  pressures,  both  with  the  foil  at  a  fixed  angle  of 
attack  and  with  it  sinusoidally  oscillating.  Analog  tape  recordings  were 
made  of  the  cavitation  noise,  angle  of  attack,  and  camera  pulse.  Figure  1 
is  a  functional  block  diagram  of  the  data  acquisition  system. 

The  cavitation  noise  was  sensed  by  a  flush  mounted  hydrophone  in 
the  top  of  the  water  tunnel's  closed  jet  test  section  46  inches  (117 
cm)  downstream  from  the  foil's  axis.  The  hydrophone's  output  was  amplified, 
filtered  through  a  four  pole  high  pass  Butterworth  filter  with  3db 
signal  attenuation  at  a  frequency  of  10  kHz,  and  then  recorded  at 
15  inches  per  second  (38cm/s)  to  IRIG  direct  record  standards.  The 
hydrophone  was  uncalibrated,  so  all  noise  measurements  were  relative 
to  an  arbitrary  level.  The  angle  of  attack,  or  foil  angle  signal,  was 
obtained  by  a  Kaman  displacement  sensor  activated  by  the  linear  deflection 
in  a  mechanical  linkage  connected  to  the  foil.  Foil  angle  calibrations 
gave  an  absolute  foil  angle  measurement  error  of  +  0.06  degrees.  The 
camera  pulse,  a  signal  generated  by  the  photographic  instrumentation, 
was  produced  each  time  a  photograph  was  taken.  The  foil  angle  signal 
and  camera  pulse  were  recorded  using  IRIG  standard  intermediate  band 

frequency  modulation  techniques.  Details  of  the  experimental  apparatus 
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and  test  procedure  are  described  by  Shen  and  Peterson. 
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127  data  runs  were  made:  32  static  runs  where  the  foil  angle  was 
fixed  at  a  given  angle  of  attack,  and  95  dynamic  runs  where  the  foil 
was  sinusoidally  oscillated  at  various  amplitudes  and  frequencies 
about  a  preset  positive  3.3  degree  attack  angle.  During  these  tests, 
the  water  speed  ranged  from  32  to  54  ft/sec  (9.7  to  16.3  m/s)  and 
the  tunnel  pressure  ranged  from  8.5  to  40  psia  (58.6  to  276  kPa) . 

For  the  static  runs,  the  foil  angle  ranged  from  a  positive  1  degree 
angle  to  positive  4  degrees.  For  the  dynamic  runs  the  foil  angle 
sine  wave  single  amplitude  varied  from  0.35  to  2.5  degrees;  a  large 
portion  of  the  runs  were  made  with  a  1  degree  oscillation.  The  oscil¬ 
lation  frequencies  were  4,  5.5,  7.5,  10,  15,  20  and  25  Hz.  The  runs 
were  divided  into  different  series;  the  tunnel  pressure  and  water 
speed  were  fixed  for  each  series.  For  the  static  runs  the  series  consisted 
of  a  range  of  foil  angle  settings.  Dynamic  series  had  the  foil  angle 
sine  wave  single  amplitude  held  constant  and  the  oscillation  frequency 
was  varied. 

HYDROPHONE  SIGNAL  CHARACTERISTICS 


The  recorded  hydrophone  signal  is  a  composite  of  background  noise, 
instrumentation  noise,  and  cavitation  noise.  Signal  processing  procedures 
must  differentiate  between  these  noise  sources  so  that  only  cavitation 
noise  is  quantified. 

General  background  noise  which  contaminates  the  cavitation  noise 
consists  mostly  of  acoustic  emissions  from  mechanical  systems  and  fluid 
flow.  The  character  of  these  emissions  changed  with  run  conditions. 
Instrumentation  noise  was  essentially  electrical  noise  from  the  tape 
recorder  • 

The  hydrophone  signal  was  high-pass  filtered  to  reduce  the  influence 
of  flow  noise.  Low-pass  filtering  was  used  to  reduce  the  influence  of  the 
tape  recorder  noise.  During  the  dynamic  runs,  the  mechanical  assembly 
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generated  noise.  Determining  the  significance  of  this  noise  source 
relative  to  cavitation  noise  required  foil  oscillation  runs  under 
noncavitating  conditions  as  well  as  identification  of  the  time  relation¬ 
ship  between  the  noise  components  and  the  angle  of  attack.  Flow  noise  was 
generated  by  turbulence  and  small  gas  bubbles  as  the  water  flowed 
past  the  foil  and  other  boundaries  in  the  test  section.  Again,  tests 
under  noncavitating  conditions,  showing  relationships  between  noise 
and  angle  of  attack  were  required  to  estimate  the  influence  of  this 
flow  noise  . 

Cavitation  noise  sensed  by  the  hydrophone  consisted  of  a  low-level 
noise  with  intense  short  duration  noise  bursts  of  random  amplitudes 
that  occurred  at  random  time  intervals.  During  dynamic  runs,  these  short 
duration  bursts  were  modulated  or  varied  by  the  foil's  changing  angle 
of  attack. 

The  variability  of  the  short  duration  noise  bursts  and  modulation 
with  angle  of  attack  were  the  two  signal  characteristics  that  prompted 
the  development  of  an  ensemble  averaging  procedure.  For  the  dynamic  runs 
ensemble  averaging  produced  the  average  sound  pressure  variance  as  a  func¬ 
tion  of  time  from  the  start  of  the  foil  angle  cycle,  i.e.  a  function  of 
foil  angle. 

DATA  REDUCTION 

This  section  gives  a  detailed  description  of  the  signal  processing 
procedures  used  to  reduce  the  experimental  oscillating  foil  data.  First, 
a  general  description  of  these  procedures  is  presented;  then  the  imple¬ 
mentation  of  these  procedures  is  detailed  in  two  parts:  Analog  Signal 
Processing  and  Digital  Signal  Processing. 


GENERAL  SIGNAL  PROCESSING  PROCEDURES 


For  each  experimental  run  the  primary  information  desired  from  the 
hydrophone  signal  was  sound  pressure  variance  (standard  deviation 
squared)  over  a  given  frequency  interval.  This  basic  measurement  took 
several  forms:  average  sound  pressure  variance  over  the  entire  run, 
sound  pressure  variance  at  each  camera  pulse  occurrence,  and  sound 
pressure  variance  as  a  function  of  foil  angle.  Since  the  hydrophone 
was  uncalibrated,  these  measurements  were  made  relative  to  an  arbitrary 
level . 

Since  cavitation  noise  varies  in  intensity  with  a  time-varying  angle 

of  attack,  the  hydrophone  signal  from  the  dynamic  runs  is  a  nonstationary 

random  variable.  Random  variables  were  treated  using  the  framework  given 

3 

by  Bendat  and  Piersol.  Each  foil  angle  cycle  was  treated  as  if  it  was  a 
new  run  at  the  same  test  conditons.  Sound  pressure  variance  was  obtained 
from  ensemble  averages  —  averages  made  from  the  square  of  hydrophone 
signal  amplitudes  measured  at  the  same  time  or  position  in  each  successive 
foil  angle  cycle.  See  Figure  2.  For  the  dynamic  runs  this  procedure 
produced  relative  sound  pressure  variance  measurements  as  a  function 
of  foil  angle. 

The  hydrophone  signal  passed  through  several  processing  phases  before 
these  relative  sound  pressure  variance  measurements  were  obtained.  First, 
the  signal  was  processed  through  an  analog  bandpass  filter.  The  choice  of 
filter  settings  qualified  the  sound  pressure  variance  measurements.  The 
next  processing  phase  produced  the  quasi-stationary  relative  sound  pressure 
variance  (QRSPV).  The  QRSPV  assumed  that  the  hydrophone  signal  was 
statistically  stationary  over  a  time  interval  of  at  least  0.0015  seconds 
for  both  static  and  dynamic  runs.  An  analog  multiplier  and  an  averaging 
filter  performed  the  operation  that  produced  the  QRSPV.  The  averaging 
filter  closely  approximated  a  running-averager^!  it  did  not  have  the  long 
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"'CUT"  AND  ORDER 
EACH  CYCLE 


QRSFvff )  =  -j-J  h2[r\dr 


f-T  T  =  0.0015  SE 


alfl^oCo+apSiNj^f 


RSPv(a|f))=  — VoRSFv(|/r-l)T 


ENSEMBLE 

AVERAGED 


"memory"  associated  with  the  standard  resistor /capacitor  (RC)  averaging 
filter.  The  exact  expression  for  the  running-averager  that  approximated 
QRSPV  as  a  function  of  time  was: 

t 

QRSPV(f)  =  rjr  jfl\)  dr 

M 


was 


where  the  integration  period  T  was  0.0015  seconds  and  f)  [  f 
the  filtered  hydrophone  signal. 

All  relative  sound  pressure  variance  measurements  were  derived  from 
this  QRSPV  time  history.  The  average  relative  sound  pressure  variance 
(ARSPV)  was  the  QRSPV  averaged  over  the  duration  of  a  nominal  30-second 
run.  The  relative  sound  pressure  variance  at  each  camera  pulse  occurrence 
(RSPVC)  was  the  QRSPV  sampled  at  each  camera  pulse.  The  relative  sound 
pressure  variance  as  a  function  of  foil  angle  (RSPV  )  was  the  QRSPV 
ensemble-averaged  over  many  foil  angle  cycles. 

Except  for  the  QRSPV  at  each  camera  pulse,  averaging  (either  run  or 
ensemble  averaging)  was  performed  on  the  QRSPV  to  build  up  statistically 
reliable  results  —  reduce  random  error  in  the  relative  sound  pressure 
variance  measurements.  The  0.0015  second  averaging  period  of  the  QRSPV 
did  not  influence  the  run  average;  however,  this  period  did  limit  or 
bias  the  detail  possible  from  the  ensemble-averaged  pressure  variance 
measurements  (trade-off  with  sample  rate  for  the  digital  processing). 

In  these  tests,  the  greatest  foil  angle  frequency  was  25  Hz,  so  this 
0.0015  second  averaging  period  effectively  divided  the  foil  angle 
cycle  into  about  27  independent  intervals. 
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The  QRSPV  signal  is  generally  quite  random  because  of  the  short 
averaging  period.  The  variability  of  this  signal  is  represented  by 
the  standard  deviation  of  the  QRSPV  divided  by  the  average  QRSPV. 

This  variability  was  calculated  for  both  run  and  ensemble  averages. 

Since  its  value  depends  on  the  averaging  filter's  integration  period, 
its  usefulness  is  restricted. 

It  does,  however,  give  a  feel  for  the  stability  of  the  noise  emission 
at  certain  foil  angles.  Also,  the  QRSPV  at  a  camera  pulse  could  be 
considered  representative  of  the  noise  usually  occurring  at  that  foil 

angle  if  its  level  is  within  +  1/2  standard  deviation  about  the  mean 
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(either  run  or  ensemble  averages).  From  random  signal  theory  the 
standard  deviation  to  mean  ratio  (random  error)  should  be  inversely 
proportional  to  the  square  root  of  the  run  length  for  gaussian  random 
processes.  Since  this  ratio  is  near  one  for  many  runs,  extrapolating 
from  the  0.0015  second  averaging  period  to  the  30-second  nominal  run 
length  gives  an  expected  random  error  below  5  percent  for  either  run 
or  ensemble-averaged  QRSPV  measurements.  These  data  are  not  necessarily 
gaussian,  so  the  extrapolation  procedure  was  checked  for  several  runs 
and  found  reasonable  for  setting  this  5  percent  error  bound.  The  check 
consisted  of  dividing  a  run  into  segments  with  lengths  1200  times 
the  QRSPV  averaging  period,  or  1.0  second.  The  QRSPV  mean  was  found 
for  each  segment.  The  standard  deviation  of  these  mean  values  divided 
by  the  mean  of  all  these  mean  values  gives  the  random  error  estimate 
for  the  1.8  second  averaging  time.  The  prescribed  random  error-to-run 
length  relationship  was  observed  by  comparing  proportionality  constants. 

Quantified  sound  pressure  variance  levels  from  statistical  run  and 
ensemble  averages  were  produced  from  a  series  of  operations  performed  on 
the  recorded  hydrophone,  foil  angle  and  camera  pulse  signals.  These 
operations  consisted  of  two  different  processes:  analog  signal  processing 
and  digital  signal  processing.  The  analog  procedures  produced  the 
QRSPV  signal,  an  amplified  foil  angle  signal  that  compensated  for  the 
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QRSPV's  running-average  phase  shift,  a  foil  angle  synchronization 
pulse  that  identified  each  positive  going  mean  level  crossing  of  the 
foil  sinusoid,  and  a  conditioned  camera  pulse.  These  data  were  digitized 
at  a  1600  Hz  sample  rate  and  then  processed  on  a  digital  computer 
to  produce  sound  pressure  variance  from  the  run  and  ensemble  averages 
as  well  as  the  sampled  values  at  each  camera  pulse.  The  details  of 
these  processing  procedures  are  presented  in  the  remainder  of  this 
report . 

ANALOG  SIGNAL  PROCESSING 


The  analog  signal  processing  procedures  produced  the  quasi¬ 
stationary  relative  sound  pressure  variance  (QRSPV) ,  phase  compensated 
foil  angle  signal,  foil  angle  synchronization  pulse,  and  conditioned 
camera  pulse.  A  block  diagram  of  the  analog  signal  processing  is  given 
in  Figure  3.  Intermediate  results  from  analog  processing  were  digitized 
and  further  processed  on  a  digital  computer  to  produce  the  final  results. 

The  recorded  data  signals  were  played  back  at  a  tape  speed  of  1-7/8 
inches  per  second  (4.8  cm/s),  1/8  their  recorded  speed.  The  lower  signal 
frequencies  simplified  the  design  of  the  analog  processing  circuitry. 
Except  where  noted,  all  references  to  time  and  frequency  made  in  this 
report  are  "real  time"  —  as  if  the  data  were  played  back  and  processed 
at  their  recorded  speed. 

QUASI -STATIONARY  RELATIVE  SOUND  PRESSURE  VARIANCE  (QRSPV) 

The  QRSPV  time  history  was  the  result  of  the  following  operations 
on  the  recorded  hydrophone  signal:  bandpass  filtering,  amplification, 
squaring  and  running  average.  The  square  root  of  the  QRSPV  was  then 
taken,  to  overcome  digitizer  limitations,  and  digitized.  The  status 
of  these  operations  was  monitored  by  a  saturation  indicator  and  a 
30-second  running  average  displayed  on  a  digital  voltmeter  (DVM). 


11 


Bandpass  filtering  was  the  first  operation  performed  on  the  hydro¬ 
phone  signal;  two  filter  configurations  were  used.  For  most  runs  this 
signal  was  filtered  by  an  Ithaco  Model  4111  four  pole  bandpass  filter. 

The  high  pass  filter  section  had  a  -3  db  attenuation  at  a  frequency  of 
8  kHz  which  reduced  flow  noise  in  the  hydrophone  signal  •  The  low  pass 
section  had  -3  db  attenuation  at  a  40  kHz  frequency  to  reject  tape 
recorder  noise.  The  effective  bandwidth  of  this  filter  was  30  kHz 
with  a  center  frequency  of  25  kHz. 

For  selected  runs  the  hydrophone  signal  was  filtered  by  a  General 
Radio  Model  1564-A  one-third  octave  filter.  The  following  center 

frequencies  were  used;  10,  12.8,  16,  20,  25.2,  32,  and  40  kHz.  From 

5 

reference  ,  the  American  Standard  Specification  for  Octave,  Half 

Octave  and  Third  Octave  Band  Filter  Sets  [ANSI  SI. 11-1966  (R  1975)], 

the  noise  bandwidth  for  these  filters  is  0.2316  f  where  f  is  the  mid- 

m  m 

band  or  center  frequency  of  the  third  octave  filter  used.  The  10  kHz 
high  pass  filter  used  during  data  acquisition  distorted  the  results 
from  the  10  kHz  center  frequency  third  octave  filter.  Measurements 
of  this  composite  transfer  function  showed  that  the  noise  bandwidth 
was  still  within  ANSI  specifications;  however,  the  QRSPV  was  reduced 
by  a  factor  of  2.06  and  its  center  frequency  was  shifted  to  10.3  kHz. 

After  filtering,  the  hydrophone  signal  was  amplified  and  then  squared. 
A  saturation  detector  monitored  the  amplifier's  output;  overrange  condi¬ 
tions  (amplifier's  output  greater  than  8  volts)  were  indicated  by  a  light 
emitting  diode  (LED).  This  indicator  was  held  on  for  about  0.5  seconds 
by  a  timing  circuit.  The  amplifier's  gain  was  selected  during  playback 
to  prevent  overrange  conditions.  These  gains  were  continuously  selectable 
from  1  to  40. 

The  amplified  signal  was  then  squared  using  a  Teledyne  Philbrick 
Model  4454  analog  multiplier.  The  square  root  function  after  the 
running-average  circuit  used  the  same  type  of  multiplier.  To  obtain 
the  squaring  and  square  root  functions,  these  multipliers  were  configured 
as  Bhown  in  the  application  section  of  the  Teledyne  Philbrick  specifi¬ 
cation  sheet  titled:  4454/4455  and  4456/4457  Accurate,  Fast  Multiply- 
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Divide-Square,  Square  Root  Operators  (15M  Revised  8/78).  The  recommended 
trim  resistors  and  set  up  procedures  were  used  with  these  multipliers. 

The  running-average  operators  used  to  obtain  the  QRSPV  (0.0015- 
second  averaging  period)  and  the  DVM  display  (30  second  averaging  period) 
give  a  continuous  output  that  closely  approximates  the  average  multiplier 
signal  present  during  the  last  averaging  period.  Signals  present  before 
that  period  have  little  influence  over  the  output.  This  averager  character¬ 
istic  tracked  the  variability  in  the  hydrophone  signal.  Examples  of  this 
performance  feature  are  shown  by  the  response  of  this  circuit  to  a  sine 
wave  burst  and  to  a  cavitation  noise  burst;  See  Figure  4.  This  performance 
was  adequate.  The  choice  of  averaging  period  was  a  tradeoff  between 
statistical  variability  in  the  QRSPV  and  the  rate  at  which  the  noise 

intensity  changed.  The  circuit  diagram  for  this  averager  (used  for 

4 

data  at  1/8  recorded  tape  speed)  is  given  in  Figure  5.  The  step  response 
of  this  circuit  (averager  alone)  is  given  in  Figure  6.  For  comparison, 
this  figure  also  contains  the  step  response  of  a  RC  averager.  The  RC 
averager  (resistor  charging  a  capacitor)  would  spread  the  effect  of  a 
large  noise  burst  over  several  time  constants  since  it  exponentially 
forgets  past  inputs. 

The  square  root  of  the  QRSPV  was  digitized.  Since  the  digitizer  had 
limited  dynamic  range  (4096  discrete  levels  cover  -10  to  10  volts),  the 
square  root  operation  was  necessary  to  compress  the  dynamic  range  of  the 
QRSPV  to  match  the  range  of  the  digitizer.  The  digitized  QRSPV  -square  - 
root  was  then  squared  during  the  digital  processing  phase  to  obtain  better 
small  signal  detail.  A  calibration  using  sinusoidal  and  random  test 
inputs  to  the  muliplier,  running-average  and  square  root  circuit  combina¬ 
tions  produced  the  QRSPV  error  curve  given  in  Figure  7.  This  error  curve 
was  used  to  estimate  the  error  in  the  final  computer  results.  The  average 
estimated  error  for  125  runs  was  7  percent  of  the  reading  for  the  run 
averaged  QRSPV.  The  highest  estimated  error  was  22  percent  for  one  of 
these  runs.  The  ensemble-averaged  QRSPV  was  generally  more  accurate  in 
the  regions  of  peak  noise  and  less  accurate  in  low  noise  regions. 
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Figure  4  -  Running  Average  Performance 
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Figure  5  -  Running  Averager  (QRSPy  &  Foil  Angle) 
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Figure  7  -  QRSPV  Error 


The  30-second  average  QRSPV  displayed  on  the  DVM  wa~  recorded  and 
later  compared  with  the  run  averaged  QRSPV  on  the  computer  output. 

These  readings  gave  preliminary  results  as  well  as  a  check  on  the 
computer  readings.  The  two  readings  were  close  even  though  they  were 
generally  obtained  from  different  sections  of  test  data  and  with  different 
run  times.  For  116  runs,  the  difference  between  the  two  readings  had 
a  standard  deviation  of  6.5  percent  of  reading. 

PHASE  COMPENSATED  FOIL  ANGLE  (PCFA) 

To  compensate  for  phase  shifts  between  QRSPV  and  foil  angle,  the  foil 
angle  signal  was  processed  through  a  running  average  circuit  with  the 
same  time  constant  as  the  unit  used  for  the  QRSPV  (0.0015  seconds). 

An  amplifier  with  a  gain  of  7  brought  the  +  1.414  volt  tape  recorder  full 
scale  signal  up  to  the  +  10  volt  full  scale  digitizer  range.  This  signal 
is  labeled  PCFA  on  Figure  3. 

FOIL  ANGLE  SYNCHRONIZATION  PULSE  (FASP) 

A  foil  angle  synchronization  pulse  (labeled  FASP  on  Figure  3)  was 
generated  for  every  positive  going  mean  level  crossing  of  the  sinusoidal 
foil  angle  signal.  This  pulse  was  digitized  and  used  in  the  computer 
processing  to  synchronize  the  ensemble-averaging  procedure. 

This  0.0015  second  duration  pulse  was  generated  by  a  one-shot  trig¬ 
gered  by  a  level  detector.  The  adjustable  level  detector  was  driven  by 
the  phase  compensated  foil  angle  signal.  To  reduce  the  influence  of  tape 
recorder  noise,  the  level  detector  had  0.2  volts  of  hysteresis  built  in; 
once  the  state  of  the  detector  changes,  the  input  voltage  must  swing  lower 
than  the  preset  level  by  0.2  volte  before  the  detector  reverts  back  to 
its  original  state.  The  level  at  which  the  detector  would  change  states 
was  adjusted  manually  by  using  an  oscilloscope  to  estimate  the  mean  of  the 
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foil  angle.  Both  QRSPV  and  foil  angle  were  ensemble-averaged  (phase  pre¬ 
served  independent  of  this  setting)  so  this  setting  was  not  critical. 

The  duration  of  the  one  shot  was  chosen  so  that  at  the  1600  Hz  digitizer 
sample  rate,  at  least  two  samples  of  the  digitized  foil  angle  synchroni¬ 
zation  pulse  would  be  at  full  pulse  voltage. 

CONDITIONED  CAMERA  PULSE  (CCP) 

The  camera  pulse  generated  when  a  photograph  was  taken  was  condi¬ 
tioned  and  digitized.  During  computer  processing,  the  QRSPV  and  foil 
angle  were  printed  whenever  one  of  these  pulses  was  detected. 

The  circuitry  used  to  condition  this  pulse  was  almost  identical  to 
the  foil  angle  synchronization  pulse  circuitry.  The  range  of  the  detec¬ 
tion  level  was  narrowed  to  improve  resolution  since  the  camera  pulse 
from  the  tape  recorder  was  a  damped  sinousoid  with  a  1  to  2  volt  peak 
value.  This  signal  is  labeled  CCP  on  Figure  3. 

An  event  counter  was  used  to  help  locate  the  camera  pulse  on  tape 
and  to  indicate  the  number  of  pulses  digitized. 

DIGITAL  SIGNAL  PROCESSING 


The  digital  signal  processing  procedures  used  the  digitized  QRSPV- 
square-root,  foil  angle  synchronization  pulse,  and  camera  pulse  to  cal¬ 
culate  the  average  sound  pressure  variance  over  the  run,  QRSPV  and  foil 
angle  at  each  camera  pulse,  and  sound  pressure  variance  as  a  function  of 
foil  angle  (ensemble  averaged  QRSPV).  Also  calculated  were  QRSPV  standard 
deviation  over  the  entire  run  and  as  a  function  of  foil  angle,  the  number 
of  camera  pulses  in  each  run,  and  the  fourier  series  of  the  ensemble  aver¬ 
aged  QRSPV  and  foil  angle. 

This  computer  processing  was  accomplished  in  three  phases.  First, 
the  time  histories  were  digitized  using  a  Hewlett  Packard  Model  2100S 
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minicomputer  and  stored  on  standard  computer  1/2  inch  magnetic  tape. 

These  tapes  were  then  processed  on  the  Center's  6000  series  CDC  digital 
computers . 

PROCESSING  PROGRAM  HYPHANS 

The  computer  program  HYPHANS  processed  the  magnetic  tapes  from  the 
HP  minicomputer  to  produce  the  desired  reduced  data.  These  results  were 
then  stored  on  permanent  file  for  use  with  the  format  program  HYPHFMT. 

The  results  are  obtained  from  one  pass  of  the  input  tape.  During 
the  pass  the  calibration  factors  are  applied  to  each  data  sample.  The 
QRSPV-square-root  samples  which  were  originally  digitized  are  squared  to 
obtain  QRSPV  samples.  These  QRSPV  samples  are  again  squared  to  produce 
QRSPV-squared  samples.  The  QRSPV  and  QRSPV-squared  samples,  and  the 
number  of  samples  are  accumulated  for  mean  and  standard  deviation 
estimates  over  the  run.  When  a  camera  pulse  is  detected,  its  location 
within  the  run,  the  foil  angle  and  QRSPV  are  stored.  For  the  dynamic 
runs,  ensemble-averaging  is  performed.  When  the  foil  angle  synchronization 
pulse  is  detected,  the  averaging  index  is  reset  and  the  foil  angle 
period  is  accumulated  for  mean  and  standard  deviation  of  the  foil 
angle  period.  As  each  time  step  is  processed,  the  foil  angle  and  QRSPV 
values  and  their  squares  are  accumulated  in  arrays  where  the  array 
element  number  is  the  current  averaging  index.  This  index  is  incremented 
by  1  for  every  time  step  that  does  not  contain  a  foil  angle  pulse. 

The  index  is  reset  to  1  when  a  foil  angle  pulse  is  detected.  Another 
array  is  used  to  accumulate  the  number  of  samples  accumulated  in  each 
array  element.  A  histogram  of  the  foil  angle  periods  is  later  extracted 
from  this  array. 

When  the  tape  pass  for  each  run  is  completed,  the  accumulated  values 
are  processed.  The  average  relative  sound  pressure  variance,  ARSPV,  over 
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the  entire  run  is  given  by: 


M 

ARSPV  =  QRSPv(fp) 

77=1 


The  QRSPV  standard  deviation  is: 


^QPSPV 


where  M  is  the  number  of  samples  in  the  run  andfy*}  is  the  time  of  each 
sample  in  the  run.  The  ensemble  averages  are  formed  using  the  same  basic 
equations  at  each  index  into  the  foil  angle  cycle.  First,  however,  an 
improved  estimate  of  the  foil  angle  period  is  made  and  the  accumulated 
values  are  again  averaged  so  that  one  foil  angle  cycle  will  always  consist 
of  64  equal  sections  or  elements  regardless  of  the  foil's  oscillating 
frequency. 
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Tape  recorder  tape  speed  variations,  noise  affecting  the  foil  angle 
synchronization  pulse  circuitry,  and  the  sampling  process  all  introduce 
time  base  jitter.  An  occasional  noise  burst  or  tape  recorder  dropout 
could  prematurely  reset  the  averaging  index.  Conditions  could  also  exist 
where  the  index  is  not  reset  until  the  end  of  the  next  cycle.  The  influ¬ 
ence  of  these  factors  on  the  final  results  must  be  reduced  or  identified. 
Ensemble-averaging  of  the  foil  angle  not  only  synchronizes  it  with  the 
ensemble  averaged  QRSPV  and  to  the  time  base,  but  it  also  provides  a 
check  on  the  ensemble-averaging  process.  It  should  produce  a  sinusoid 

with  a  known  amplitude  and  period  to  compare  with  previously  processed 

1 

experimental  data.  The  foil  a^gle  period  histogram  also  aids  in  verifying 
the  test  results.  A  refined  estimate  of  the  foil  angle  period  is  made 
by  taking  the  original  "raw"  values  for  the  period's  mean  and  standard 
deviation  (sigma)  and  then  averaging  this  histogram  over  the  interval 
from  mean  minus  4  sigma  to  mean  plus  4  sigma.  The  raw  and  refined 
values  are  printed  in  the  output  to  verify  the  procedure. 

These  steps  were  taken  to  improve  the  quality  of  the  data;  however, 
little  change  was  noted.  Only  about  five  of  the  runs  had  their  period 
modified  by  more  than  one  sample  period  (reciprocal  of  sample  rate).  The 
final  period  standard  deviation  (time  base  jitter)  was  generally  less 
than  1  percent  of  the  mean  foil  angle  period.  The  foil  angle  amplitudes 
and  frequencies  were  compared  to  previously  processed  results  where  power 
spectrum  methods  were  used  to  estimate  these  signal  parameters.  The  stan¬ 
dard  deviation  of  the  difference  between  the  two  amplitudes  was  the  equi¬ 
valent  of  0.01  degrees.  The  standard  deviation  of  the  difference  between 
the  two  frequencies  was  0.016  Hz  for  a  frequency  range  from  4  to  25  Hz. 

The  ensemble  averages  are  formed  by  dividing  this  refined  foil  angle 
period  into  64  increments  and  accumulating  (interpolating  as  necessary) 
the  ensemble  accumulations  over  these  new  increments.  Equations  similar 
to  those  presented  earlier  for  ARSPV  and  \y  are  then  applied  to 

these  accumulations  at  each  increment  where  M  is  also  an  accumulated  value. 
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Once  these  values  are  obtained,  they  are  printed  and  stored  on  per¬ 
manent  file.  With  this  program,  the  output  format  Is  not  very  descriptive. 
Another  program,  HYPHFMT,  presents  the  results  In  a  more  useful  format. 

FORMATTING  PROGRAM  HYPHFMT 

The  program  HYPHFMT  takes  the  results  stored  on  permanent  file  by 
HYPHANS,  combines  them  with  a  test  condition  summary  file,  and  reproduces 
them  on  output.  The  Fourier  series  of  ensemble-averaged  foil  angle  and 
QRSPV  Is  also  presented,  as  are  plots  of  ensemble-averaged  QRSPV.  The 
Appendix  is  a  sample  output  from  HYPHFMT.  Run  number,  tunnel  conditions 
(water  speed  and  tunnel  pressure),  nominal  foil  angle  and  amplitude, 
and  the  measured  foil  angle  parameters  (amplitude,  frequency  and  mean) 
from  previous  processing  are  then  combined  with  the  HYPHFMT  outputs. 

Four  pages  of  computer  processed  results  are  obtained  for  each  run . 

The  top  of  each  page  contains  the  test  conditions:  run  number  and  digital 
tape  file  number;  tunnel  water  speed  and  pressure;  nominal  foil  angle 
amplitude  and  frequency;  the  nondimensional  coefficients  of  Reynolds  number, 
cavitation  number  and  reduced  frequency;  and  the  foil  angle  single  amplitude 
and  mean  from  the  earlier  data  reduction. 

In  addition,  the  first  page  contains  the  calibration  correction 
factor,  sample  rate,  hydrophone  variance  or  QRSPV  averaged  over  the  run, 
and  the  camera  pulse  results.  The  sample  rate  presented  in  the  output 
is  the  rate  used  at  1/8  the  recorded  speed.  The  hydrophone  variance 
or  QRSPV  results  contain  the  number  of  samples  in  the  run,  the  mean  or 
average  QRSPV  over  the  run,  the  QRSPV  standard  deviation  and  the  standard 
deviatlon-to-mean  ratio.  The  camera  pulse  results  are  next:  the  element 
or  pulse  number,  the  period  in  samples  between  camera  pulses  (the  number 
for  the  first  element  is  from  start  of  run),  the  foil  angle  at  that  pulse 
and  the  hydrophone  variance  or  QRSPV  for  that  pulse. 


The  second  page  contains  tables  of  the  ensemble-averaged  (signal 
averaging)  results*  First,  the  total  number  of  foil  angle  pulses 
and  the  foil  angle  period  and  standard  deviation  are  printed  from 
both  the  "raw"  and  the  refined  (adjusted)  procedures.  The  ensemble- 
averaged  table  contains:  the  element  from  1  to  64  (foil  angle  period 
divided  into  64  intervals) ,  the  number  of  samples  used  to  produce 
the  averages,  the  ensemble-averaged  foil  angle  mean  and  standard  deviation, 
the  hydrophone  QRSPV  (variance)  mean  and  standard  deviation,  and  the  QRSPV 
standard  deviation-to-mean  ratio.  The  other  five  columns  contain  the 
discrete  Fourier  series  results:  the  harmonic  index  and  the  magnitude 
and  phase  for  ensemble-averaged  foil  angle  and  QRSPV. 

The  discrete  Fourier  series  of  the  ensemble-averaged  foil  angle  and 
QRSPV  is  generated  by  using  Fast  Fourier  Transform  (FFT)  procedures.  The 
magnitude  and  phase  terms  presented  in  the  output  are  the  terms  necessary 
to  reproduce  ensemble-averaged  foil  angle  and  QRSPV  at  each  sample  point. 
The  following  equation  is  used: 

3' 

ar)=<3o  +  L  iK  siN|ff/<1rM+0K) 
k=i  . 

where  Q^are  the  magnitude  terms  and  are  the  phase  terms  at  each 
harmonic,  K  .  Complete  reproducibility  is,  however,  not  possible  because 
the  32nd  harmonic  was  not  saved  during  the  FFT  computations. 

The  third  and  fourth  output  pages  contain  line  printer  plots  of 
ensemble-averaged  QRSPV  and  its  standard  deviation-to-mean  ratio.  For 
each  of  the  64  points  these  results  are  scaled  by  100  increments;  the 
maximum  value  generally  sets  the  value  of  the  100th  Increment.  The  scaling 
interval  is  printed  at  the  top  of  the  plot.  The  increment  value  at  each 
foil  angle  is  printed  in  the  left  column.  The  values  for  these  plots 
come  directly  from  the  tables  on  the  second  page. 


CONCLUSION 


This  application  of  ensemble-averaging  principles  is  a  good  example  of 
their  power  in  satisfying  the  data  reduction  needs  for  certain  classes 
of  variable  or  nonstationary  data.  The  results  produced  by  these 
procedures  were  used  by  Shen  and  Peterson  to  analyze  the  cavitation 
noise  generated  by  an  oscillating  hydrofoil.  These  quantitative  results 
allowed  the  correlation  of  cavitation  noise  with  observed  phenomena 
such  as  cavitation  Inception  and  termination  and  cavity  length.  The 
magnitude  and  location  within  the  foil  angle  cycle  of  the  generated 
noise  were  important  elements  in  this  analysis. 

These  data  reduction  procedures  produced  good  results.  The  run 
lengths  were  long  enough  to  keep  statistical  sampling  error  (random 
error)  less  then  5%.  The  circuitry  used  to  produce  the  QRSPV  measurement 
was  generally  operated  at  signal  voltage  levels  that  resulted  in  an 
average  error  of  7%.  And  time  base  jitter  was  generally  less  than 
IX  of  the  foil  angle  period.  Since  the  foil  angle  amplitude  and  fre¬ 
quency  obtained  from  the  ensemble  averaging  procedure  gave  values  that 
were  nearly  identical  to  earlier  results  obtained  using  power  spectrum 
methods,  the  ensemble  averaging  procedure  was  well  verified. 

Use  of  these  procedures  could  benefit  other  experimental  programs 
in  which  random  data  is  modulated  by  a  cyclic  time  function.  Noise  ex¬ 
periments  with  cavitating  propellers  is  one  such  example. 

Additional  Improvements  could  be  made.  Dedicated  hardware  and  a 
minicomputer  at  the  test  site  would  allow  near  real  time  response.  Also, 
since  the  analog  multipliers  introduce  most  of  the  error  associated  with 
these  procedures,  improvement  is  possible  with  better  equipment  selection. 
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